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The Swift Gamma-Ray Burst redshift distribution: 
selection biases and optical brightness evolution at high-z? 
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ABSTRACT 

We employ realistic constraints on astrophysical and instrumental selection effects to 
model the Gamma-Ray Burst (GRB) redshift distribution using Swift triggered red- 
shift samples acquired from optical afterglows (OA) and the TOUGH survey. Models 
for the Malmquist bias, redshift desert, and the fraction of afterglows missing because 
of host galaxy dust extinction, are used to show how the "true" GRB redshift distri- 
bution is distorted to its presently observed biased distribution. We also investigate 
another selection effect arising from a correlation between E{so and Lopt. The analy- 
sis, which accounts for the missing fraction of redshifts in the two data subsets, shows 
that a combination of selection effects (both instrumental and astrophysical) can de- 
scribe the observed GRB redshift distribution. Furthermore, the observed distribution 
is compatible with a GRB rate evolution that tracks the global SFR, although the rate 
at high- 2: cannot be constrained with confidence. Taking selection effects into account, 
it is not necessary to invoke high-energy GRB luminosity evolution with redshift to 
explain the observed GRB rate at high- 2; . 

Key words: gamma-rays: bursts - (ISM:) dust, extinction - methods: statistical 



1 INTRODUCTION 

Long duration gamma ray bursts (GRBs) provide a po- 
tentially powerful probe of the early Univer se. In particu- 
lar, t he spatial distribution of Swift satellite (|Gehrels et al.l 
I2OO4I ') GRBs could be used to infer the star formation 
rate history at 2: > 2. Although Swift has localized nearly 
800 GRBs up to late 2011, only about 50% of these have 
detected optical afterglows (OAs), despite rapid optical 
follow-up. Furthermore, only about 30% have measured red- 
shifts. Despite this deficit, there have been numerous at- 
tempts at analysing/simulating the GRB redshift distribu- 
tion to constrain the star formation rate and/or GRB lumi 



Yiiksel et all I2OO8I: IWanderman fc PiranI I2OI0I : lOin et all 



2010l : IVirgili et aLlboill ) 



Previous work using the GRB redshift distri- 
butio n as a constraint on rate/luminosity evolution 
/e.g. iKocevski fc Lianel I2OO6I : IWanderman fc PiranI I2OI0I : 
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20lj ) has not addressed the problem of 
redshifts. Given that only ~ 30% of Swift 
the spatial distribution of the miss- 



nosity evolution at h i gh-z (see e.g. Nataraian et al, 
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the "missing" 
bursts have redshifts, 
ing redshifts plays a fundamental role in understanding the 
link between GRBs, their environments, and relationship to 
the star formation rate. The first critical step in identify- 
ing where the missing redshifts are located is to understand 
why GRB optical afterglows are not observed for ~ 5 0% of 
Swift bursts. In a previous studv lc"oward et al.l (|2008l ) iden- 
tified and constrained a global GRB redshift selection effect 
function, but used a purely empirical approach that did not 
include optical imaging limits. 

In this study, we use realistic constraints and models for 
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redshift dependant selection biases, combined with GRB OA 
luminosities, to show how selection effects distort the "trure" 
spatial distribution to its presently observed distribution. 
We employ two subsets of GRB redshifts. The first, Howell 
& Coward (2012), hereafter HC, uses 141 Swift triggered 

fectroscopic absorption redshifts from OAs up to Oct 2012. 
The second, less biased but smaller sample, uses a subset 
of 58 redshifts from the TOUGH (The Optically Unbiased 
GRB Host) survey (|Hiorth et al.|[201^ ). By accounting for 
selection effects, we investigate if the observed GRB redshift 
distribution is compatible with GRB rate evolution tracking 
the global star formation rate. 



2 GRB OPTICAL SELECTION EFFECTS 

2.1 Instrumental selection effects 

We define GRB OA selection effects as the combination of 
sensitivity limited optical follow-up and phenomena (astro- 
physical and environmental) that reduce the detection prob- 
ability of an OA. Below we list the main instrumental biases 
which inhibit both OA detection an d redshift measurement. 
These are discussed bv lFvnbo et al. I (|200a ) and lHiorth et al.l 
l|2012h in more detail. 

(i) XRT localisation: An XRT localisation improves the 
probability of OA detection, especially in crowded fields. 

(ii) Foreground extinction: Galactic dust extinction 
significantly reduces the probability of detecting an OA. 

(iii) Source declination: The lack of rapid follow-up in- 
struments capable of imaging at 5 < ±70° . 

(iv) Source angular distance from Sun: Bursts that 
occur too near the Sun have a limited follow-up window in 
time. 

Historically, because of the deficiency in pre-Swift 
ground-based follow-up of GRBs, there was a strong bias for 
imaging the brightest bursts. Because the brightest bursts 
are predominantly nearby, a significant faction of the first 
GRB redshifts were obtained by emission spectroscopy of the 
host galaxy. In both the pre-Swift and Swift era (from 2005 
onwards), an optical afterglow (OA) is usually required to 
measure a redshift. For most high-jz GRBs, this is achieved 
by absorption spectroscopy of the GRB afterglow. The host 
galaxies are usually too faint to make a significant contri- 
bution to the spectra. Most GRB spectroscopic redshifts 
are acquired by ground based te lescopes, including VLT, 
Gemini-S-N, Keck and Lick (see iFvnbo et al.l |2009| . for a 
more complete list along with specific spectroscopy instru- 
ments) . 

The measurement of a GRB redshift depends strongly 
on the limiting sensitivity and spectral coverage of the spec- 
troscopic system. For example, the VLT FORS 2 instrument 
has limiting magnitude of ~ 24, while the VLT X-Shooter 
has a limit of m ~ 22 for 3 x 60 min exposure times. One of 
the challenges for GRB spectroscopy is the fact that the 
OA is fading, so the delay time for acquiring spectra is 
also a factor. This bias is expected to manifest at high-z, 
where the brightest OAs are near the limiting sensitivity 



^ GRB redshift sample is a subset taken from GCN circulars and 
|http://www.mpe.mpg.de/~jcg/grbgen.html| 
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Figure 1. The Swift mean redshift uncertainty bound plotted 
over the duration of the mission. It is clear there is a drift in the 
mean redshift over time, a consequence of different priorities and 
instruments con tribut ing to redshift acquisition i.e. the learning 
curve effect (see ICow ard 2009). The jump observed in 2009 is a 
result of GRBs 090423 and 090429B, with redshifts of ^ = 8.26 
(NIR spectroscopic) and z = 9.2 (photometric) respectively. 



for spectroscopy. For faint OAs, photometric redshifts are 
measured by multi-band opt ical/NIR telescopes. For exam- 
ple, GROND Ig reiner et al. I [2008) is providing photometric 
redshifts of high-z GRBs via a seven channel simultaneous 
imager mounted on the 2.2 m MPl/ESO telescope at the 
ESO/LaSilla observatory. 

Other selection effects are more subtle, for example im- 
proving instrumentation and observation techniques over 
time that are used for both the OA follow-up and redshift 
measurement (|Cowardll2009l ). The learning curve effect can 
be considered as a time-dependant bias because it affects 
the redshift distribution differently over time. An example 
is the slow drift in the median redshift during the first 4 
years of Swift operation. Fig. [T] shows the trend towards a 
smaller Swift mean redshift over the duration of the mission. 



2.2 Redshift dependant selection effects 

This following analysis focusses on the remaining selection 
effects that distort the redshift distribution over certain red- 
shift ranges, with the assumption the GRB OA brightness 
is a proxy for redshift measurement i.e. a reasonably bright 
OA is usually required to obtain an absorption redshift. 



(i) Malmquist bias: This bias arises because the tele- 
scopes and instruments acquiring OA absorption spectra 
(and photometry) are limited by sensitivity. In reality, the 
instruments acquiring redshifts are biased to sampling the 
bright end of the OA luminosity function. To account for this 
bias, it is necessary to have some knowledge of OA luminos- 
ity function (which is uncertain especially at the faint end), 
and an estimate of the average sensitivity limit of the instru- 
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Figure 2. The normalized redshift distribution of bright QSOs 
with mji < 20 (blue curve), faint QSOs m/j > 20 (black curve), 
and very faint QSOs mn > 20.9 (red curve) from the Sloan Digital 
Sky Survey Quasar Catalog (7th Data Release). It is clear that 
the selection of fainter QSOs results in increasing selection effects. 
It is more difficult to acquire redshifts in the redshift desert (0.8 < 
z <2) for the faintest sources. 



in z = 2 — 2.5, which are more suited for follow-up spec- 
troscopy in the near-IR. 

(iii) Host galaxy extinction: There has been growing 
evidence that dark bursts arc obscured in their host galaxies 
(e.g. [jakobsson ct al. 2004; Levan ct al. 2006; Pcrlcy ct al] 
I2OO9I : ISvensson et al.l 120121 ). The detection of the near-IR 
OAs of some GRBs (which would have been considered 
as dark bursts because their OAs were not detected in 
any bluer bands) provides evidence for dust obscuration 
(e.g. Tanvir et al. 2008). These studies generally show that 
GRBs originating in very red host galaxies always show 
some evidence of dust extinction in their afterglows. 

IPerlev et al.l (|2009l ) show that a significant fraction of 
dark burst hosts have extinction column s with Ay ~ 1 mag , 
and some as high as Ay = 2 — 6 mag. I Rossi et"al] (|2012l ) 
performed a search for the host galaxies of 17 bursts with 
no optical afterglow. They find in seven cases extremely red 
objects in the error circles, at least four of them might be 
dust-enshrouded galaxies. The most extreme case is the host 
of GRB 080207, one of the reddest galaxies ever associated 
with a GRB. Using mul t i- wave length observations of GRB 
080207, ISvensson et all (|2012l ) show that that the dark 
GRB hosts are systematically more massive than those 
hosting optically bright events, perhaps implying that 
previous host samples are severely biased by the exclusion 
of dark events. 



ments. This is the most fundamental bias that encompasses 
all peak fiux limited detection and is the basis for modelling 
a selection function for OA/redshift measurement. We also 
note another selection effect for very high- 2 GRBs. For very 
high-z GRBs the Malmquist bias is partially reduced by time 
dilation that extends the emissions implying it is easier to 
measure a redshift. 

(ii) Redshift desert: The so-called redshift desert is a 
region in redshift (1.3 < z < 3) where it is difficult to mea- 
sure absorption and emission spectra. This becomes more 
significant for faint sources (small signal to noise ratio spec- 
tra). As redshift increases beyond z ^ 1, the main rest-frame 
optical emission features from the host galaxies (as for ex- 
ample the [OUT] A3727 line, H/3, Ha) become more difficult 
to detect since they move to the IR, a wavelength region 
where the sensitivity of CCDs starts to drop and sky bright- 
ness increases. At the same time, most of the strong UV 
absorption features in the OA spectra as those produced by 
the low ionization state metals (e.g. Fell, SII, Sill), and the 
weak rest-frame UV lines of active star forming galaxies are 
not redshifted enough to be detected less than z — 2 — 2.5. 

To demonstrate the effect of source brightness on num- 
ber counts. Fig. [2] shows the redshift distribution of faint 
quasars from the Sloan Digital Sky Survey Quasar Catalog 
(7th Data Release). For the faintest sources, (i.e. mu > 20.9) 
it is more difficult to acquire spectroscopic redshifts because 
these sources are most affected by the redshift desert, as 
shown by the strong 'dip' in the redshift distribution. 

It is difficult to quantify the fraction of missed desert red- 
shifts across different surveys because o f the varying strate - 
gies and instruments used. For example, ISteidel et al.l (|2005l ) 
find a typical success rate of (70-80)% for obtaining spectro- 
scopic redshifts of photometrically selected galaxies in the 
desert. This success rate is expected to improve for galaxies 



IZafar et all (|201ll ) compared the high redshift GRBs to a 
sample of lower redshift GRB extinctions and found a lack 
of even modera t ely ex tinguished events {Av ~ 0.3) at z > 4. 
iBouwens et~aLl (|2009l ) studied the effect of dust obscuration 
on estimates for the high-z star formation rate. Their work 
showed that star-forming galaxies at z > 5 almost univer- 
sally have very blue UV-continuum slopes, and that there 
are not likely to be a substantial number of dust-obscured 
galaxies at 2 > 5 that are missed in "dropout" searches. 
In spite of the biased selection and small number statis- 
tics, their analysis supports the argument for a decrease in 
dust content in star-forming environments at high redshifts. 
Hence it is highly likely that dust obscuration in GRB host 
galaxies is linked to the evolving star formation rate. This 
implies that the missing redshifts (with obscured OAs) may 
be tracking the star formation rate of red dusty star forming 
galaxies, while the observed redshift distribution is tracking 
the star formation rate of bluer more compact star forming 
galaxies. Supporting the hig h-dust at small redshift scenario, 
IWatson fc JakobssonI (|2012l ) show that GRBs with high X- 
ray absorption column densities found at 2 < 4 typically 
have very high dust extinction column densities, while those 
found at the highest redshifts do not. 

The most significant issue with interpreting a lack of 
highly extinguished OAs at high-2 as evidence for less dust 
extinction, is selection effects. Specifically, the non-obscured, 
bright OAs are preferentially selected at high-2, so that a 
possible obscured subset of GRBs remains hidden. Our in- 
dependent approach to constrain an upper limit on GRB 
optical dropouts from host galaxy dust with redshift is 
based on optical sur veys of core-collapse supernovae (SNe). 
iMattila et al.l l|2012l ) find for a volume-limited rest-frame op- 
tical SN survey, the missing SN fraction increases from its 
average local value of about 19% to 38% at z — 1.2, and then 



© 3002 RAS, MNRAS OOO.mfTOl 



4 



1 1 1 1 1 1 1 1 1 1 1 1 1 — 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 





Figure 3. A statis tical comparison of thre e sub-sets of GRB re d- 
shifts: from Salvat erra et al'l ll20lj') - (S), iHiorth et al.l hOlt i ~ 
TOUGH and the subset of lHowell fc CowardI ||2013| ') "ThO. plot- 
ted in Figure [4] The HC sample has 3 times as many bursts, and 
is a very similar distribution to that of TOUGH; Pks=94%. Al- 
ternatively the S sample, with a peak flux cut-off of 2.6 ph s~l, 
is biased against high-2 redshifts. 

decreases marginally beyond z = 2. iMattila et al.l (|2012l l 
note that their estimated missing SNe fraction is also con - 
sistent with the recent findings of iMichalowski et al.l (|2012l ^ 
for their sample of 2 < 1 GRB host galaxies. 



2.3 Data selection 

Our primary redshift sample is based on the Howell and 
Coward (2012) subset from the J. Greiner catalogue of well- 
localized GRBs consisting of 141 Swift triggered spectro- 
scopic absorption redshifts from OAs up to Oct 2012. We 
note that the HC redshifts are mostly measured from the 
OA, as opposed to the host galaxy, so that OA brightness 
is the key factor in determining a redshift. We also include 
photometric redshifts with 2 > 6, because of the difficulty 
for obtaining spectroscopic redshifts at very high-2 because 
of Lyman-Qf dropout of optical sources. 

We also compare two other catalogues of high redshift 
sub-samples o f Swift long GRBs: the TOUGH sample and 
the sample of ISalvaterra et al] l|2012f ). which we label as S. 
For comparison we show in Figure [3] the HC sample (141) 
with the samples of TOUGH (53 bursts) and S (51 bursts): 
we note that in this latter sample 4 bursts (GRB050802, 
GRB060306, GRB060614, GRB081221) have uncertain red- 
shifts according to the HC sample and an additional burst, 
050416A, is catagorised as an XRF. Figure shows that the 
the HC sample has a very similar distribution to that of 
TOUGH, i.e. a two sample Pks=94%. The S sample shows 
that a peak flux cutoff at 2.6 ph s~l has produced a signif- 
icant deficiency of bursts above 2 ~ 2 in comparison with 
the TOUGH and HC samples. This is supported by the cor- 
responding two sample Pks probabilities of 9% and 2% re- 
spectively. We find that by applying similar peak flux cut in 
our sample (reducing our sample to 52 bursts) our distribu- 
tion resembles that of S, yielding a Pks probability of 83% 
(and 45% when compared to the TOUGH sample). 

We choose to employ the HC subset because it: 
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Figure 4. Top: The redshift distribution of 147 Swift GRBs 
(up Oct 2012) acquired by absorption spectroscopy of the optical 
afterglow, the HC sample, (normalized bars), but including six 

2 > 6 NIR spectroscopic and photometric redshifts, along with 
a number density curve (dashed line) that smoothes small scale 
fluctuations. Bottom: The TOUGH survey subset of 58 redshifts 
obtained using strict selection criteria to minimize observational 
biases. 

a) demonstrates high statistical compatibility with the 
TOUGH sample; b) allows us to use a larger sample (nearly 

3 fold number increase compared to the TOUGH and S 
samples). We conclude that the technique of removing the 
faintest GRB peak fluxes, sample S, introduces a redshift 
dependant bias that manifests as a deficiency of high-2 red- 
shifts, as shown in Figure [31 

The HC redshift subset is representative of about 37% 
of the total Swift sample. The TOUGH survey subset of 
58 redshifts, obtained using strict selection criteria to mini- 
mize observational biases, is 55% complete in redshift. This 
means that the two data sets require different normaliza- 
tions to account for the missing redshift fractions. For the 
TOUGH data set, we assume that the instrumental biases 
are mostly accounted for, so that the remaining 45% of red- 
shifts are missing because of the three redshift dependant 
biases, namely the redshift desert, host galaxy dust obscu- 
ration and the Malmquist bias. For the HC subset, the nor- 
malization is 18% missed because of instrumental biases, and 
55% missed because of the redshift desert, host galaxy dust 
obscuration and the Malmquist bias. Figure [4] plots the HC 
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and TOUGH redshift distribution along with their density 
approximations. 



3 GRB REDSHIFT DISTRIBUTION MODEL 
WITH SELECTION EFFECTS 

The dominant GRB redshift distribution biases discussed 
above are represented as the product of independent dimen- 
sionless selection functions that are unity for a 100% selec- 
tion probability: 

(i) V'Obs - number dropouts from mostly non-redshift de- 
pendant biases listed in i\2.1\ which are different depending 
on the selection criteria for the sample. We assume that the 
TOUGH sample is relatively free of instrumental biases, but 
about 20% of the HC sample is affected by instrumental bi- 
ases. 

(ii) V'Swift(z) - the limited sensitivity of Swift to trigger 
on GRBs (see Appendix [XI for a derivation of tlisvUtiz)) ■ 

(iii) tpM{z) - the limited sensitivity of instruments to mea- 
sure a redshift from the GRB OA. 

(iv) VDcsort(2) - number dropouts from the redshift 
desert. 

(v) 7/)Dust(z) - number dropouts from host galaxy dust 
extinction. 

The GRB redshift probability distribution function, 
that includes the above selection effects, can be expressed 
as: 



P(z) = Np^:^^-^^i>Sv,iit{z)tpObsipM{z) 



dz (1 + z) 
cscrt (z)^Dust(2:) 



(1) 



where is a normalization constant. The volume element, 
dV/dz, is calculated using a flat-A cosmology with Ho = 
71 km s~^ Mpc~^, Q.M = 0.3 and Q.h ~ 0.7, and we fix 
V'Obs ~ 0.5 (see the selection effects listed in H2.ip . The 
function e[z) is the dimensionless source rate density evo- 
lution function (scaled so that e(0) = 1). We assume that 
e{z) tracks the star formation rate history, and we employ 
the star formati on rate best fit to multi- wavelength surveys 
from Fig. 10 in Red dv &: Steidell (|2009ll , shown in Fig. [S] 
Because the fit only extends to 2 = 6, we extrapolate the fit 
to z = 10. Fig. lBll plots e{z) using a simple functional form: 
two linear increases from z = Q — 2.5, and an exponential 
decrease a.t z > 2.5 (see Appendix IB)) . 

The question of how the GRB rate evolves w ith metal- 
licity remains controversial. lEUiott et all (|2012l ) found no 
evidence for GRB rate evolution tracking metallicity from 
a sub-sample of GR OND, and that various redshift biases 
l|Coward et al.l [20081 ) could account for previous interpreta- 
tions. Because of this unceratinty, we adopt this argument 
and do not include metallicity evolution, but instead assume 
GRBs track the global evolving star formation rate, i.e. e{z) 
above. 

The total selection function, i.e. excluding ipsv,m{z), is 
defined as: 

^Total = -i/'Obs X iIjm{z) X l/>Dcsort(2) X ?/>Dust (z) . (2) 

The total fraction of GRB redshifts that are missed because 
of optical selection effects can be expressed as 

(■10 



N-'CP{z)dz 



Sl° dV{z)/dz{l + z)-ie(2)^swift(2)dz 




Figure 5. Top panel: Star formation rate history from multi- 
wavelength surveys taken from Fig. 10 in lReddv &: Steidell ||2009| ) 
and reproduced by permission of the AAS. We use the solid 
line, for the best-fit star-formation history assuming a luminosity- 
dependent dust correction to z 2. See Appendix [B] for the con- 
version between the SFR and GRB rate evolution model, e{z). 



where the normalization is cancelled to convert to a rela- 
tive number density. This provides an independent consis- 
tency check of V'Totai , because Ft ~ 63% (45%) for HC and 
TOUGH respectively, is based on the fraction of measured 
redshifts to the total number of Swift bursts. We apply the 
Kolmogorov-Smirnov test (KS test) test using the model 
P{z) that includes selection effects, and the two data sets 
HC and TOUGH. In addition to the KS test we employ an 
error calculation of the fraction of redshifts missed in the 
model compared to the observed missed fraction. The frac- 
tional error, E is: 



E= {Ft + I - N)/N, 



(4) 



where I is the fraction missed because of instrumental bi- 
ases, which we assume is zero for TOUGH and 0.18 for the 
HC sample. The normalization constant, A^, for HC and 
TOUGH is 0.63 and 0.45 respectively. 



3.1 Optical afterglow selection function 

The GRB OA luminosity function (LF) in R-band is ap- 
proximated by fi tting to the compiled OA luminosities from 
iKann et al.l l|2010 ). Their luminosities, scaled to 1 day, have 
been corrected for Galactic extinction, host extinction and 
are k-corrected if possible. We use the functio nal form of 
IJohannesson. Biornsson fc GudmundssonI l|2007l ) to approx- 
imate the LF: 



,(L) = C(|- 



exp 



\n^L/Lo) 
2cr2 



exp 



(3) 



where C is a normalization constant, Lo a characteristic lu- 
minosity. We convert L/Lq to absolute magnitudes using 
L/Lo = io(Afo-M)/2.5 ugjjjg ^Yie KS test as a constraint 
{Pks > 0.8), we approximate the LF using Mo — 24.5, 
a = 3.5 and A = 0.5. This LF was converted into a prob- 
ability distribution function (PDF), plotted in Fig. [T] by 
normalizing over the luminosity range [—28, —16], where we 
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have extended the faint end to M = —16 to crudely account 
for the dropout of faint OAs from the Malmquist bias. 

We note that the above LF does not take into account 
the Lyman-a cut-off at high-z R-band sources. Lyman-a 
absorption renders sources at z > 6 undetectable in the R- 
band. Nonetheless, the GRB redshift data contains 6 sources 
with either spectroscopic or photometric redshifts at 2 > 6, 
which are significant in our analysis. For these very high-z 
GRBs, we assume that the OA brightness in filters redder 
than R-band (including the near-infrared) follows approxi- 
mately the same number distribution as in R-band (ahhough 
in the H-filter the GRB OAs will be brighter from 1-2 mags 
but the CCD sensitivity reduces by a similar magnitude). 

3.1.1 A Malmqmst bias tmth a Eiso — iopt correlation 

In a study of bo th short and long duration GRBs 
iGehrels et al.l l|2008t ) found optical, X-ray and 7-ray emis- 
sions to be linearly correlated. A highly significant correla- 
tion (> 99.99%) was found between the prompt 7-ray fiuence 
and X-ray flux of long duration GRBs (r 0.53 ). A slightly 
less significant correlation (99%) was revealed between the 
optical and X-ray fluxes (correlation coefficient, r ^ 0.44). 
Similar correlations existed for the short duration GRB sam- 
ple, although were less significant due to the smaller sample. 
Stronger bursts in the X-ray and 7-ray produced brighter af- 
terglows in the optical and X-ray respectively. A large spread 
was seen around the c orrelation line in all instances. 

iKann et al.l (|2010l ) flnd a weak correlation between iJiso 
and Lopt with a Kendall rank correlation coefficient of 0.29 
for the subset of GRBs used in this study to estimate an op- 
tical LF. The key question is: how does this correlation affect 
the Malmquist bias? Firstly consider the effect if Eiso and 
^opt are uncorrelated. This implies that a high _Eiso (that 
will be preferentially detected by Swift) could be associated 
with equal probability with either a low or high luminosity 
optical afterglow. This scenario implies that detection of the 
OA is independent of the high energy luminosity. Alterna- 
tively, observation suggests a high Eiso preferentially selects 
a high optical luminosity. This implies that the Malmquist 
bias will be reduced at high-z, because it is the high Eiso 
bursts that are preferentially seen at large-z. These bursts 
will also be more optically luminous, so that redshift mea- 
surement will be more probable. 

Figure E| plot s the 1-day optical luminosities (at z=l) 
(|Kann et alTfeoiol 'l (79 bursts) against bolometric isotropic 
burst energy. We also plot the sub-sample used in this study 
(from the HC sample) - obtaining a subset of 48 bursts. On 
the same Figure we calculate the change in mean optical 
luminosity across -Eiso by binning the data. It is clear that 
there is a trend for increasing mean OA luminosity with 
increasing -Eiso. The change in magnitude, AM, across is 
-Eiso is about AM « 1.9. As described above, the correlation 
will reduce the Malmquist bias at high-2:. We employ the 
constraint of a maximum increase in luminosity of AA/ ~ 
1.4 mag in z = — 6. The Mamquist correction is defined as 
AM{z) = 1.7Logi(,(l + z), so that AM(6) fa 1.4. 

3.1.2 Malmquist bias 

A telescope limiting magnitude (mi,) defines a threshold for 
obtaining a redshift. Because of the different instruments 
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Figure 6. Top panel: The 1-day luminosities (at z=l) for 
the Kann sample of optical afterglows (79 bursts) are shown as 
symbols against bolometric isotropic burst energy. The sub- 
sample used in this study - 48 bursts obtained via spectroscopic 
absorption and designated as LGRBs in the HC catalogue are in- 
dicated by diamonds. We find a weak correlation with a Kendall 
rank correlation coefficient of 0.21 (with 97% significance) - this 
is less than 0.29 obtained by Kann et al (2010) for their full op- 
tical luminosity sample. Bottom panel: The mean magnitudes 
and their l-c errors obtained per logarithmically spaced bins of 
equal width. We use 7 bins (indicated by horizontal the error 
bars) which is the maximum number that obeys a requirement of 
at least 5 events per bin. The drift in the mean magnitude with 
increasing -Eiso of the sample, AM, is about 1.4 mag. 
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Figure 7. The GRB OA probability distribution function (PDF) 
using equation \5\ and GRB OA luminosities from Kann et al. 
(2010). The optimal model is constrained by the KS test (Pks > 
0.8) with Mq = 24.5, cr = 3.5 and A = 0.5, and normalized over 
the Magnitude range [—28,-16]. We note that the faint end is 
the least constrained because of the Malmquist bias. 
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engaged in spectroscopic redshift measurement it is dilBcult 
to constrain this important parameter. For definiteness, we 
take TTiL ~ 24, approximating an average limit for the tele- 
scopes used to acquire the redshifts for the faintest bursts. 
We include an average extinction from Kann et al. (2010) of 
Av = 0.3. 

With a reference time tc of 1 day and assuming that all 
OAs fade in luminosity by , the OA limiting luminosity 



was calculated at a time Tz (|lmerito et al.|[2009l ) as: 

MlW =mL-51ogio (^)-|logio (^'^yk{z)+AM{z)~Av. 



For t he mean time taken to acquire a redshift JCowardl 
l2009t ). we use an average Tz ~ 400 min. Hence Tz links the 
probability of redshift measurement to the brightness of the 
OA. For the k- correction w e use — 2.5(/j — l)Login [1 + ^], with 
/? — 0.66 from lKann et al.1 (|2010| ). We employ two scenarios 
for the limiting magnitude. The first assumes a correction for 
the correlation of higher Eiso with brighter optical luminosi- 
ties, AM{z) (see section[6l). We also calculate the Malmquist 
bias without this correction. Using the above definitions, the 
OA selection function defines the fraction of observable OAs 
with redshift: 

iPoa{z) = / ip{M)dz . (7) 

The total fraction of GRB redshifts missed because of 
the Malmquist bias is less than 10% using the above model 
assumptions. 

3.2 OA host galaxy dust obscuration 

For an estimate of optical GRB numb er dropouts from host 
galaxy dust, we employ the study of iMattila et aD (j2012l l 
using a volume-limited rest-frame core-collapse SN survey. 
Their missing SNe fraction is obtained in two parts: First, 
they estimate the fraction of SNe in normal galaxies with 
substantially higher host galaxy extinctions than predicted 
by simple models for the smooth dust distribution and the 
resulting inclination effects. 

Secondly, they estimate the fraction of SNe missed in 
local U/LIRGs. For this they make use of the rich SN popu- 
lation of one of the nearest LIRGs, Arp 299. Assuming that 
the SNe with the highest host galaxy extinctions are missed 
by the optical searches and not compensated for by the stan- 
dard extinction corrections, they derive the fraction of SNe 
that remain missing and estimate the corrections needed to 
be applied when deriving CCSN rates. Finally they use this 
information together with the latest knowledge of the nature 
and evolution of high-z U/LIRGs to calculate the fraction 
of SNc misse d as a funct ion of redshift. 

[Mattila et al.l (|2012h provide a parameterization of the 
form k X z + m for the missing SNe fraction as a function 
of redshift out to z = 2. We extend this parameterization 
beyond z = 2 by extrapolating the function beyond z = 2. 
Using this model for the fraction of missed GRB OAs from 
host galaxy dust extinction, we obtain a missing GRB OA 
fraction of 30-35%, which is comparable to recent estimates 
based on spectroscopy of GRB host galaxies. We note that 
this model assumes that GRB OAs are either detected (be- 
cause of minimal host extinction) or not observed at all be- 
cause of very high host extinction. 



0.5 




0.5 



■A M correction 
■ no A M correction 



4 6 
redsliift 



10 



Figure 8. Top panel: The red shift selection function, 'i/'Dust(z) 
is based on lMattila et al. I ll2012h . who provide a parameterization 
for the missing SNe fraction as a function of redshift out to z = 2. 
We extend this parameterization beyond z = 2 hy extrapolating 
the function values from 2 = 2 — 10. Middle panel: The curve 
shows the selection function for the redshift desert, based on a 
85% success rate of redshift measurement in galaxy surveys in 
2 = 1.3 — 2.5. Bottom panel: The redshift selection function 
iPm{z) for two scenarios, the first, with rriL = 24, a mean time 
for redshift measurement Tz of 400 min, and a correction term 
for the Eiso — ^optical correlation. The second scenario, dashed 
curve, is the same as above but without this correction. 



3.3 Redshift desert 

In addition to the effect of dust extinction, the redshift 
desert decreases the probability of acquiring a spectroscopic 
redshift in the 1.3 < z < 2.5 interval. Quantifying this prob- 
ability within this region is very difficult, because it depends 
on the strategy and brightness of the sources. As an obser- 
vational constraint on the probability of redshift measure- 
ment in the desert we employ the succes s rate of a cquir - 
ing galaxy spectroscopic redshifts from Ste idel et al.l (|2005l l 
i.e. (70 — 80)%. The simplest functional form of i/'Dcsort (-z) 
is a smoothed step function with tpDcscit{z) — 80% for 
1.3 < z < 2.5 and unity probability outside this region. 
Fig. [8] plots the selection effect models ipoustiz), V'Desert(z) 
and tl^Miz). The total fraction missed over 2 = — 10 is less 
than 10%. 



4 RESULTS AND DISCUSSION 

The inclusion of selection effects is required to achieve the 
maximum K-S probability, for the two data sets (HC and 
TOUGH), using the redshift cumulative distribution and 
model. Fig. [9] plots the observed redshift distribution, with 
the optimal model (Model 1 below), and for comparison the 
expected distribution that would be observed if all optical 
selection effects were removed. Below we list the four model 
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Table 1. The maximum KS statistics (Pks) using the model cu- 
mulative distributions for the HC and TOUGH samples. Models 
with selection effects (Models 1, 2 and 3) assume mL = 24 and 
= 400 min. Also shown for each Model are the total fraction 
of unmeasured redshifts, Fj-, and the error calculated from the 
Swift fraction of redshifts. 





Pks 


e(2 = lO)"!" 


AA/(z)t 




_Fy error% 


HC 












Model 1 


0.86 


1 


yes 


65% 


3% 


Model 2 


0.27 


4 


yes 


71% 


13% 


Model 3 


0.96 


1 


no 


70% 


11% 


Model 4 


0.13 


1 


no 


0% 


100% 


TOUGH 












Model 1 


0.68 


1 


yes 


47% 


4% 


Model 2 


0.17 


4 


yes 


53% 


18% 


Model 3 


0.86 


1 


no 


52% 


16% 


Model 4 


0.15 


1 


no 


0% 


100% 



t GRB rate evolution e{z 
the fit of lReddv 
AM(z) ^ 1.7Logio(l + z) 



10) = 1 is an extrapolation from 
Steidell ll2009h . ^ the Malmquist correction is 



scenarios and Table 1 tabulates the K-S probabilities and 
fractional errors for the models: 

(i) Model 1 - Includes selection effects and the 
Malmquist bias correction term {/S.M{z)), assuming jtil = 
24 and e{z — 10) — 1, i.e. a SFR at z = 10 similar to that 
at z = 0. 

(ii) Model 2 - Includes selection effects, excludes 
Malmquist bias correction term, with e{z = 10) = 4, i.e. 
a SFR at 2 = 10 four times higher than at z = 0. 

(iii) Model 3 - Includes selection effects, excludes 
Malmquist bias correction term and uses the same GRB 
rate evolution as Model 1, e{z — 10) = 1. 

(iv) Model 4 - Excludes all optical selection effects. This 
model does not account for the missing redshifts from the 
HC and TOUGH samples. 

Even though model 4 (which excludes all optical selec- 
tion effects) is not rejected by the KS test {P = 13%), it 
does not account in any way for the missing 45% and 63% 
(TOUGH and HC data) of redshifts from the Swift sample. 
Alternatively, Model 1, including realistic selection effects, is 
compatible with the observed distribution including a GRB 
rate evolution that follows the global SFR. 

It is useful to stress how different selection effects could 
have similar affects on the redshift distribution. For example, 
there are different biases that could potentially reproduce 
excess number counts at high-z: 

(i) A larger intrinsic GRB rate at very high-2 {z — 10) 
via SFR or metallicity evolution compared to the SFR at 
z = (Model 2). 

(ii) Less dust obscuration of GRB hosts at high- 2;, result- 
ing in brighter OAs. 

(iii) GRB circumburst environments could be relatively 
different at high-z, which would also affect the OA bright- 
ness. 

(iv) Correlation between high Eiso and Lopt bursts boost 
the fraction of redshifts measurable at high-2. 

(v) An evolving GRB high energy LF. Given the correla- 



tion between Eiso and Lopt, an evolving high energy lumi- 
nosity will further increase the optical luminosities at high-z, 
so that redshift measurement is more probable. 

In summary, our analysis suggests that a combination of 
selection effects (both instrumental and astrophysical) can 
adequately describe the observed redshift distribution. Fur- 
thermore the observed distribution is compatible with a rate 
evolution that tracks the evolving SFR. We have also high- 
lighted that two GRB redshift catalogues claiming to be 
statistically complete, are in fact not compatible with each 
other. We show that the TOUGH selection and a subset of 
absorption redshifts (the HC sub-sample) are compatible. 

The optimal strategy to avoid the complicated selec- 
tion biases in the current redshift sample is to employ dedi- 
cated ground based instruments with well understood sensi- 
tivities and performance. Robotic instruments with integral 
field low resolution spectroscopic capabilities could also play 
an important role, both in terms of constraining redshifts 
and studying the early optical emissions. Also, high-energy 
detection combined with simultaneous optical coverage will 
continue to play a role in increasing the number of redshifts, 
especially when combined with a dedicated ground-based 
NIR follow-up network. 



ACKNOWLEDGMENTS 

D.M. Coward is supported by an Australian Research Coun- 
cil Future Fellowship. E.J. Howell acknowledges support 
from a University of Western Australia Fellowship. 



REFERENCES 

Band D.L., 2003, ApJ, 588, 945 
Bouwens R. J., et al., 2009, ApJ, 705, 936 
Butler N. R., Bloom J. S., Poznanski D., 2010, ApJ , 711, 
495 

Coward D.M., Guetta G., Burman R.R., Imerito A., MN- 

RAS, 2008, 386, 111 
Coward D.M., 2009, MNRAS, 393, L65 
Daigne F., Rossi E. M., Mochkovitch R., 2006, MNRAS, 

372, 1034 

Elliott J., Greiner J., Khochfar S., Schady P., Johnson J. L., 

Rau A., 2012, A&A , 539, A113 
Fynbo J. P. U., et al., 2009, ApJS, 185, 526 
Gehrels N., Chincarini C, Giommi P. et al. 2004, ApJ, 611, 

1005 

Gehrels N., et al., 2008, ApJ, 689, 1161 
Greiner J., et al., 2008, PASP, 120, 405 
Hjorth J., et al., 2012, ApJ, 756, 187 
Howell, E. J., & Coward, D. M. 2013, MNRAS , 428, 167 
Imerito A., Coward D., Burman R., Blair D., 2008, MN- 
RAS, 391, 405 

Imerito A., Coward D.M., Burman R.R., Blair D.G., 2009, 

MNRAS, 399, L108 
Kocevski D., West A. A., Modjaz M., 2009, ApJ, 702, 377 
Jakobsson P., et al. 2004, A&A, 427, 785 
Jakobsson P., et al., 2012, ApJ, 752, 62 
Johannesson G., Bjornsson G., Gudmundsson D.H., 2007, 

A&A, 472, L29 
Kann D. A., et al., 2010, ApJ, 720, 1513 



© 3002 RAS, MNRAS OOO.nifTOl 



9 



with selection effects 
no selection effects 
GRB data 




10 



with selection effects 
no selection effects 
GRB data 




4 6 
redshift 



10 
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model, Model 4, excludes selection effects. Bottom: Same as the 
top figure but using the TOUGH redshift distribution. Both the 
HC and TOUGH data require the same optimal models (1 and 
3). 
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APPENDIX A: 
FUNCTION 



THE SWIFT GRB SELECTION 



The GRB selection function is defined as the fraction of po- 
tentially observable bursts by Swift as a function of redshift. 
It depends on the LF of the bursts and the sensitivity of the 
BAT detector to GRBs. To account for the flux limit of the 
Swift BAT detector we include two corrections: the cosmo- 
logical k-correction, and Bdet, the bolometric correction of 
the detector: 



rlOOOOkcV 
JlkcV 



EN{E)AE 



n ENiE)dE 



(Al) 



where the interval [ei , 62] is the spectral energy band that 
the detector is sensitive to i.e (150— 15)keV and N(E) is the 
ensemble photon flux density (phs~^ cm~^ keV~^) at the 
photon energy E. The model b urst spectrum we employ is 
the Band function (lBandll2003l '). with a = -1 and /3 = -2.3 
and a rest frame break energy of 511 keV. We employ a 
k-correction of the form: 



k{z) 



I'S^'' ENiE)dE 



I< 



(l+z)150kcV 
{l-|-z)15keV 



EN{E)dE 



(A2) 



For the GRB high energy LF we use th e single power 
law form used by IPorciani fc Madaul l|2000t ) which has an 
exponential cutoff at low luminosity: 

■i'{L) = <fo \^v[-^) ■ (A3) 

Here, L is the peak rest frame photon luminosity in the 1- 
10000 keV energy range, a ensures an asymptotic slope at 
the bright end and L, is a characteris t ic cut off scaling. This 
form was used by iHowell fc CowardI l|2013l l to perform x 
minimisation of the Swift LGRBs brightness distribution. 
We adopt their best fitting parameters of: L, = 2.01q q2 ^ 
10^^ erg s"^ and a = 3.8115:6. 

The fraction of detectable GRBs, or flux-limited selec- 
tion function (those observed with peak flux > Fum) is 

V'Swift(^) = / 0(i)dL, (A4) 
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Figure Bl. Assuming tliat GRB rate evoluti o n trac ks tlie SFR, 
we extrapolate the fit from iReddv fc Steidell ||2009|'| to ^ = 10, 
using equation IIBll l. We convert this fit to a dimensionless evo- 
lution function, e{z) =SFRD{2)/SFRD{0). The solid curve is the 
extrapolation to 2; = 10 used in Model 1, and the dashed curve 
is used to produce Model 2 (see Appendix [B] for the parameters 
used for e{z)). 



where the luminosity is an isotropic equivalent luminosity 
and -Liim(-Fiim, z) is obtained by solving 

Jlim(^lim,2) = , ^iT, . k{z)Bdet, (A5) 

where di^{z) denotes the luminosity distance of the burst 
using a flat-A cosmology with Ho = 71 km s~^ Mpc~^, Qm 
= 0.3 and flA = 0.7. The flux limit, Fum, for GRB detection 
of Swift, we use the smalles t peak flux of 0.04 ph cm~^ s~^ 
from the lHowell fc CowardI (|2013l ) sample. 



APPENDIX B: GRB RATE EVOLUTION 

Assuming GRB rate evolution tracks the SFR, we use the 
following function to approximate GRB rate evolution: 

Uz z < 0.8 

e(z) = I 7 + 22 0.8 s; 2 < 2.5 (Bl) 

a[exp(2 - S)]"*" 2>2.5 

For the optimal GRB redshift distribution model plotted in 
Fig. [9] we employ a = 9.5 and b = 0.42. For the model ex- 
cluding OA brightness evolution, the optimal model requires 
e{z) parameters a = 11 and b — 0.25, yielding a GRB rate 
evolution about 4 times higher at 2 = 10 compared to the 
present day (see Fig. IB1|I . 
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